Astronomy & Astrophysics manuscript no. 
(will be inserted by hand later) 



On Shadow of the Moon in Extensive Air 

Shower Data 

F. Sheidaei'''*, A. Anvari"'^'**, M. Bahmanabadi"''''***, and M. Khakian 

Ghomi'^'^^ J. Samimi"''''-'- 

a) ALBORZ Observatory, Sharif University of Technology, Tehran, Iran. 

b) Department of physics, Sharif University of Technology, P.O.Box 11365- 
9161, Tehran, Iran. 

the date of receipt and acceptance should be inserted later 

Abstract. A new technique has been devised for the analysis of extensive air 
shower data in observing the effect of the moon on this data. In this technique 
the number of EAS events with arrival directions falling in error circles centered 
about the moving moon is compared to the mean number of events falling in 
error circles with centers randomly chosen in the sky. For any assumed angular 
radius of the error circle the deficit in EAS event count in the direction of moon 
which is a moon-related effect is interpreted as the shadow of the moon. A simple 
theoretical model has been developed to relate Nsky to the angular radius of the 
error circle and has been applied to the counts from the moon's direction in order 
to extract the physical parameters of the shadow of the moon. The technique and 
the theoretical model has been used on 1.7 x 10^ EAS events recorded at Alborz 
observatory. 
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1. Introduction 

In 1957, Clark (jClark 1975^ recognized that the moon and the sun cast a shadow in the 
isotropic flux of cosmic-ray nuclei. As they pass overhead during a transit, the moon 
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and the sun block cosmic rays, so their shadows in the cosmic ray flux should be vis- 
ible to Extensive Air Shower (EAS) arrays with sufflciently good angular resolution. 
It was not until the early 1990's that an extensive air shower array had the required 
angular resolution to observe this shadow. Since that time the shadow of the sun and 
the moon have been used to measure the angular resolution of extensive air shower 
arrays (|Urban et al. 1990| [AiEbrosio et al. 20031 lAmenomori et al. 2007p . Obse rving the 
shadow of the moon in EAS experiments which usually might have a much larger error 
circle than the disk of the moon is a very difficult task and requires a careful scrutiniza- 
tion of the data. Normally study on the shadow of moon is needed as the first step to 
find angular resolution and systematic errors of any observational data which is to be 
used for any astronomical study in particular for the research of point sources in sky. 
The object of this paper is to describe a new technique and procedure devised for the 
analysis of extensive air shower data in observing the effect of the moon on this data 
and also to present a theoretical model from which we can extract the relevant physical 
parameters of the shadow of moon. The technique has been previously briefly outlined 
on an investigation of energy spectrum of EGRET 7— ray point sources with EAS ex- 
periment (jKhakian et al. 2005 ap . In Sec. 2 we will describe this new technique and its 
procedure for EAS data analysis. In Sec. 3 we will present a theoretical model from which 
we can extract the relevant physical parameters, In Sec. 4 the result of applying the proce- 
dure and the theoretical model to 1.7 x 10^ EAS data of Alborz observatory is presented 
and the obtained results on shadow of the moon is also presented. Sec. 5 is devoted to 
discussion and some concluding remarks. 

2. Description of the New Technique and EAS data Analysis procedure 

This technique and the corresponding EAS data analysis procedure for observing the 
effect of the moon on the data and determining the pertinent physical parameters is 
based on corrected measured EAS data (corrected for systematic errors). The data must 
contain the following information for each EAS event: The arrival time of each shower, 
tg, and the coordinates of the arrival direction of the shower. The local coordinates 
of each shower arrival direction should be converted to declination and right ascension 
of each shower with arrived direction denoted here by Ss , RA^ respectively. Here, we 
denote the declination and right ascension of the center of the moving moon by 6m (t) and 
RAm{t) at the arrival time of each shower. The EAS experiments generally might have an 
experimental uncertainty which results in the experiment's error circle to be larger than 
the disk of the moon. Furthermore, normally, realistic determination of the radius of the 
error circle is best obtained by observation of the shadow moon in EAS data. However, 
in this new proposed technique and the procedure described here, there is no need for 
the use of a predetermined radius of error circle or a fitted value for it, and, instead it is 
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based on the observation of any possible deficit in shower counts falhng in the error circle 
centered about the moving location of moon as compared to the average shower counts 
falling in error circles centered at other positions in sky during the observation time for 
a wide range of assumed values for the radius of the error circle ranging from 0.2° to a 
maximum value relevant to the particular EAS data set under analysis. For any assumed 
radius of circle of error, or equivalently, its angular radius 9err the number of showers 
falling in each circle is determined by calculating the angular separations 9sm between 
the arrival direction of each shower event {5s , RAg ) and the direction of the center of the 
moon at the time of recording of that event, using the following equation from spherical 
geometry: 

cos 6sm = cos 6m COS Sg + sin 6m sin 6s cos{RAm — RAg). (1) 

Obviously, if 9sm < 9err that shower is counted as falling in the moon's error circle. In 
order to compare the obtained result with random sampling and scrutinize the difference 
for each assumed value of Oerr some random locations in the sky denoted by celestial 
coordinates {6r,RAr) arc chosen and the number of showers falling in the error circles 
centered about each of the random locations is determined similarly by calculating the 
angular separation Ogr of each shower arrival direction {Ss,RAs) with the direction of 
the center of the randomly chosen error circle denoted by {6r, RAr) from above equation 
with {6m, RAm) replaced by {6r, RAr). If for any shower event Osr < Oerr that shower is 
counted as falling in the error circle of that random position. For any assumed angular 
error radius, Oerr, some error circles are chosen in the sky centered about truly random 
locations. The number of random circles for every angular error radius should be limited 
such that it ensures that no overlap occurs between two or more random circles, so the 
number of random centers are varied from at least 1000 (for small error circles) to 77 a 
smaller number which depends on the data set under analysis (for larger error circles). 
Thus, for each assumed radius for the error circle, the mean of the shower counts falling 
in random circles could safely be used as the expected mean number of EAS events 
falling in the error circle in any direction in the moonless sky against which the number 
of EAS falling in the circle centered about the moving moon could be safely compared, 
and, the variance of its distribution could safely be used as the statistical error of the 
mean number. Obviously if the deficit in the number of showers falling in each error 
circle centered about the moving moon from the mean number exceeds the statistical 
uncertainly in the mean , then the deficit could only be attributed to the moon's effect. 
This effect is the shadow of moon in the EAS data, and, as explained in Sec. 3, from a 
quantitative analysis and comparison with the expected mean number of EAS showers 
falling in the randomly centered error circle with that falling in the circles centered about 
the moving moon, the physical parameters of the moon's shadow could be extracted. 
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3. A theoretical model for Moon's shadow in EAS data 

Following the procedure described in above technique, one now has the expected mean 
number of showers falling in randomly chosen error circles in monless sky as well as that 
in the moon-centered error circle for a set of assumed radii of the error circles, Perr-, 
(or equivalently angular radii, 6 err)- Here we propose a simple physical model to obtain 
this expected mean number as a function of the angular radius of the error circle, as 
a function of 6 err- We now derive the expected mean number of EAS events from each 
random direction in the moonless sky as a function of the assumed radius for the error 
circle, Nsky{Perr)- The derivation is based on the single assumption of the model, that 
is, the assumption of uniform intensity, I, of EAS producing radiation everywhere in the 
47r steradian for each area of the random moonless sky and thus for each element of the 
error circle. For the contribution from each element of area of the error circle (27rpdp), 
we should take into account only a fraction of radiations coming at such an angle to be 
able to reach the point of observation, that is a fraction equal to Where dQ. is the 
solid angle subtended by the element of area element from the observation point, which is 
the projection of the area element divided by the square of its distance from observation 
point, that is, dO. = 2'Kpdp{^)^, with R = (? -\- and d is the distance from point 
of observation to the center of the error circle, and it is merely a multiplicative constant 
factor relating the radius of the error circle to its angular radius Perr = d tan ^err- The 
integration is trivial. Thus we have: 

for the number of showers falling in the error circle centered about moon, Nmooniferr) 
the integration has to be split into two or tree parts involving the physical parameters 
of shadow of moon in EAS data. Here, we define the following three physical parameters 
used in this model: 

a) rm = radius of umbra of shadow that is, from p ^ to p = rm the EAS producing 
radiation are assumed to be fully absorbed (totally blocked) and could not contribute to 

A^ moon^Perr ) • 

b) Tp = radius of penumbra of shadow; that is from p = to r = Vp only a fraction (/) 
of EAS producing radiation penetrate the penumbra and contribute to N^nooniPerr)- 

c) / = the fraction of EAS producing radiation which penetrate the moon's penumbra. 
Obviously, if perr < fp the integration will only be split into two parts, that is, Jq"'^^ 

ox/;'"+./x/;-. 

For the case of perr > rp, the integral will be split into three parts: 

^ X +/ X +1 X . 

The result of these integrations giving NmooniPerr) in terms of physical parameters of 
the moon's shadow {rm,rp,,f) is given in Appendix. It should be emphasized here that 
the strict explicit assumption of uniform flux of EAS producing radiation used in this 
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model requires that when the EAS data is used to extract the parameters of the moon's 
shadow from this model one has to make sure that the data may only have statistical 
error (as this has also been implicitly assumed as outlined in the procedure for obtaining 
shower counts in the error circles), that is, the data should have been corrected for any 
systematic errors such as those related to the site of observation and non- Uniformities 
in the exposure time in various directions in sky. 



4. Application of the Technique to ALBORZ EAS data 

The technique described in Sec. 2 for determination of the moon shadow has been ap- 
plied to 1.7 X 10^ EAS data collected in 280 hours of observations in April- June 2002, 
with the small EAS array of the prototype of ALBORZ Observatory of Sharif University 
located in Tehran, Iran (51° 20'e and 35° 43'n, elevation 1200 m = 890 g cm^^). For 
details of array and data, see (jBahmanabadi et al. 2003^ . As explained in our previous 
report (jKhakian 2005 bp . the data has been corrected (scaled for uniform exposure) for 
site-dependent factors effecting shower counts from different directions in sky. The in- 
formation on the celestial coordinates of the moon during the observation time of the 



collected data has been obtained from the internet site (http://aa.usno.navy.mil). The 
moon's data has been obtained for time increments of one minute, and the location of 
moon in Right ascension and declination coordinate at the recorded time of arrival of 
each EAS event has been calculated and used. 

For an assumed set of radii for error circles ranging from 0° to 6° (increments of 0.2°) 
we have calculated the number of corrected shower counts falling in each error circle. For 
every assumed radius, some random moon-like locations passing through the paths like 
as moon's path in sky was chosen according to the radius of error circles and the number 
of corrected EAS counts falling in each error circle was found according to the procedure 
described in Sec. 2. In table 1 the number of random circles and the mean of counts for 
various assumed radii of error circles are given in second and third columns, also in this 
table (4th column) the number of corrected shower counts from the error circle centered 
about the moving moon is shown. In column 5 and 6 of the table the deficit of counts 
from moon's error circle from the mean count of the random sky error circle and the 
statistical error of deficits are given. The last column of the table gives the statistical 
significance of these deficits calculated with Li&Ma method IjLifcMa 1983P . Fig. 1 shows 
the variation of the mean number of events for moon- like error circles with random cen- 
ters as a function of the chosen radius of the error circles. The smooth curve shown that 
calculated according to our theoretical model of Sec. 3 and it fits the mean count from 
random sky with a regression of 0.996. Error bars are taken from 4th column of table 1. 
The good fit of random sky counts with model shows that we can safely use these mean 
number of events to compare with that falling in the error circles centered about the 
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moving moon and rule out the possibility of the deficit in the number of events falling 
in the moon centered circle as due to statistical fluctuations. In Fig. 2 we have shown the 
variation of events falling in each circle centered about the moon and the mean number 
of events falling in the error circles centered about random moving moon-like locations 
as a function of radius of the error circles. As seen in Fig. 2 moon counts are less than 
mean counts from random moon-like centers for all error circles that we considered. In 
Fig. 3 we have shown the number of deficit events for each radius of the error circles. We 
have fitted the deficit counts falling in the error circle centered about the moving moon 
from that for moon-like circles with random centers to our theoretical model (sec. 3 and 
Eq. A2 in Appendix) and have obtained the following results: 
0,„ = O.5°,0p = 4.5°,/ = 0.80. 

5. Concluding Remarks 

It is worth remarking that the application of the proposed technique to ALBORZ 
EAS data has yielded good agreement between the mean number of counts from error 
circles with centers chosen randomly in sky with no moon in the line of sight and the 
expected number according to our proposed theoretical model. This good agreement 
is very encouraging and prompted us to extract the physical parameters of the moon 
shadow (defined in Sec. 3) from this data. It should also be remarked that the data used 
for calculating shower counts in each error circle was the corrected counts scaled in 
order to obtain a uniform exposure of sky. The correction accounted for site-dependent 
systematic errors arising from uneven number of EAS events in various directions in sky 
due to two main factors: (1) varying amount of air mass which produces the EAS event 
as a function of zenith angle and depends on the elevation of the site (jKhakian 2005 b|) . 
and (2) geomagnetic effect which depends on the components of magnetic field at the 
site's location. Our attempt to extract the physical parameters of moon's shadow from 
this data has been fully successful as can be seen from the reported result in Sec. 4. That 
is in fitting the corrected data to our theoretical model we are able to obtain a value for 
the radius of shadow's umbra 9m = 0.5°. However , according to statistical significance 
shown in the last column of table 1 we didn't see the umbra with good significance 
but the obtained results show that in spite of low-statistics EAS data base this method 
is powerful to find shadow of moon. One may suggest that the value of 9p = 4.5° we 
obtained is just the umbra's radius rather than penumbra's and resulted from low 
angular resolution of our array . This could be right since the angular resolution of 
our array which was reported before (jKhakian et al. 2005 a|) is about 4.3° close to 4.5° 
which we find here as the radius of penumbra. We believe that the main uncertainty in 
extracting results from Alborz data could be due to the following two reasons, both of 
which will be improved upon in future with much higher number of events and larger 
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statistics. 

1. Inaccurate and low-statistics EAS data base 

Since the umbra's radius is in the order of 0.5° it is hard to expect to extract it from 
inaccurate EAS data. The local coordinates associated with each EAS event in ALBORZ 
EAS data has been obtained from an array with a very small number of detectors. 
EAS data from observatories with large arrays, once corrected for the systematic 
site-dependent errors may be more accurate to yield better results for the physical 
parameters of moon's shadow according to the technique presented here. 

2. Incomplete data on Moon 

As explained in Sec. 3 the information about the celestial coordinates of the moon 
was obtained from the internet site using time increments of one minute. The time of 
arrival of EAS events had been recorded with an uncertainty of 0.07 seconds. In the 
computations of shower counts in the error circles of various radii centered about the 
moving moon which arc given in Table 1 to check whether a given EAS event falls in the 
error circle centered about the moon or falls outside it. we have used the coordinates of 
the moon at the one of the minute steps which is closest to the arrival time of the given 
EAS event. Obviously, this may have caused an extra inaccuracy in the counts given in 
column 5 of Table 1. In future application of this technique the interpolated or exact 
location of the moon at the instant of recording of each EAS event must be used and 
the variable earth-moon distance should also be taken into account. However, the study 
of the moon's motion has been beyond the scope of the present work. 



6. Appendix 

For calciilating the count in the error circle centered about the moon, we split Eqn.2 in 
three parts. The result of integrations for two regions arc shown in following equations: 

Nmoon = -{[{C0s9err " COsOp) + f{cOs{9p - C0s6'„i))], 9err > Op 

Nmoon = —|-/(C0S (?err — C0s9„i), 9^ < 9err < 9p (^-1) 

If we now subtract the number of events in the random circles, Ngi^y from above 

we obtain: 

^back ^moon — 2 [cOS ^e^.^. 1] ^err ^ 

Nback - Nmoon = " /) C0s9err - 1 + fcOs9m] 9rn < 9err < Op {A.2) 

Nback - Nmoon = -| [(1 - /) COS 6*^ - 1 -|- / COS 9m] 9p < 9err 

The parameter /(=497270) is determined by fitting Eqn.2 to the data of column 3 in 
table 1. By knowing / and fitting the data of deficit events (column 6 of table 1) with 
above equations (A.2), we obtained O^^Op, and/. 
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Table 1. Number of EAS events obtained in various error circles with random centers 



and moon center, of the low-statistics EAS data of Alborz observatory. 
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Fig. 1. Variation of the mean number of EAS events falling in the random error circles 
as a function of the angular radius of the error circle, 9err- The smooth curve is the result 
of computations according to our theoretical model (Sec. 2). 
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Fig. 2. Variation of events falling in the random circles(») and moving moon(A) as a 
function of the angular radius of the error circle, the smooth curve is the same as Fig.l. 
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Fig. 3. Variation of deficit events falling in moving moon circles from that in the random 
circles as a function of the angular radius of the error circle, the smooth two parts curve 
is the result of fitting data with the theoretical model equations see in appendix. 



